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Abstract 

We probe the SU(3) vacuum using eigenvectors of the Dirac op- 
erator with an arbitrary phase for the temporal boundary condition. 
We consider configurations with topological charge \Q\ = 1 near the 
QCD phase transition and at low temperatures on a torus. For all our 
ensembles we show that the zero-mode of the Dirac operator changes 
its position as one changes the phase of the boundary condition. For 
ensembles near the QCD phase transition our results closely resemble 
the behavior of zero-modes for Kraan - van Baal solutions of the clas- 
sical Yang-Mills equations where the individual lumps are interpreted 
as monopoles. Our findings near T c and on the torus show that for 
both cases an excitation with topological charge \Q\ = 1 is built from 
several separate lumps. 
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1 Introduction 



QCD exhibits two remarkable features, confinement and chiral symmetry 
breaking, which both are cornerstones for our understanding of strongly 
interacting particles. It is an interesting property of QCD that at the QCD 
phase transition the restoration of chiral symmetry and the disappearance of 
confinement occur at the same critical temperature. This seems to indicate 
that the two phenomena are intimately linked with each other. So far the 
two phenomena have resisted our attempts to find a description unifying 
them and the structure of the fundamental excitations of the QCD vacuum 
is still unclear. 

In recent years lattice calculations have started to contribute to our 
understanding of the QCD vacuum fluctuations. In particular it was found 
that the eigenmodes of the Dirac operator provide a powerful filter removing 
the hard UV fluctuations and only the long range modes couple to the low 
eigenvectors These studies were inspired by the phenomenological 

picture of chiral symmetry breaking based on instantons jHj ■ In this picture 
the QCD vacuum is described as a fluid of interacting instantons and anti- 
instantons. Instead of a zero-mode which exists for a single instanton Ej 
the fluid of interacting instantons and anti-instantons gives rise to many near 
zero-modes which have small but non-zero eigenvalues. The density of these 
eigenvalues is then related to the chiral condensate by the Banks-Casher 
formula |Hj. Lattice studies have established that the near zero- modes do 
indeed show a lumpy structure and are locally chiral as expected from the 
instanton picture. Furthermore it was demonstrated that the lumps in the 
field strength have a high degree of self-duality 3 . Note however, that first 
attempts to fit eigenmodes of the Dirac operator with the profile of the 't 
Hooft zero-mode have failed 

In this article we use a new method for probing the QCD vacuum with 
eigenvectors of the lattice Dirac operator. We introduce an arbitrary phase 
exp(i2-7r£) with £ G [0, 1] at the temporal boundary of the Dirac operator 
and analyze how the zero-mode in configurations with topological charge 
Q = ±1 responds to changes in £. The most prominent feature we discover 
is that the zero-mode changes its position and can be located at different 
space-time points for different values of £. This property is observed for 
quenched ensembles of SU(3) gauge configurations with temperature, both 
in the confined and the deconfined phase as well as for configurations on a 
torus. 

Our generalized boundary condition is motivated by an interesting prop- 
erty of zero-modes for Kraan-van Baal (KvB) solutions of the classical Yang- 
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Mills equations on a euclidean cylinder. KvB solutions generalize the 
caloron solution ^U] by allowing for non-trivial Polyakov loop at spatial in- 
finity KvB solutions depend in addition to the phases of the Polyakov loop 
at spatial infinity also on N (for SU(iV)) vectors yi,i = 1,2... N. When 
one draws apart those vectors one finds that an object of charge 1 is built 
from N constituent monopoles and this property might lead to the missing 
link between confinement and chiral symmetry breaking. Strong evidence 
for SU(2) KvB solutions in cooled lattice gauge configurations were given 
for twisted jllj and periodic boundary conditions |12j . 

For KvB solutions also the zero-mode has been computed [T^]- This 
zero-mode has the remarkable property that it is located on only one of 
the constituent monopoles but it can jump from one monopole to another 
when changing Q. In a previous article ^3] we have compared zero- modes 
with periodic and anti-periodic boundary conditions and have found clear 
signals for KvB-type behavior. Here we now use the generalized boundary 
condition with phase exp(i27r£) and present further strong evidence that the 
excitations of QCD at high temperature have the structure of KvB solutions. 

In a subsequent step we apply the techniques that were successful in 
detecting the constituent monopoles at high temperature to gauge configu- 
rations generated on a torus. For the torus no equivalent of KvB solutions is 
known. However, several articles have put forward the idea that also at low 
temperatures constituents with fractional charge build up lumps with inte- 
ger topological charge [T^]- Our finding, that when changing the fermionic 
boundary condition the zero-mode does change its position, indicates that 
indeed an excitation with topological charge \Q\ = 1 is built from several 
separate lumps. 

2 Technicalities 

2.1 Gauge configurations 

For our runs we use gauge ensembles generated with the Liischer-Weisz 
action ^B] with coefficients from tadpole improved perturbation theory. We 
work at two different values of the inverse gauge coupling (3 = 8.20 and 
(3 = 8.45. A determination of the lattice spacing based on the Sommer 
parameter gives a = 0.115(1) fm and a = 0.094(1) respectively |17j . 

At those two couplings we generate configurations on 6 x 20 3 lattices. 
The two values of (3 then give rise to an ensemble in the confining phase 
(/? = 8.20) and an ensemble in the deconfined phase {(3 = 8.45). The 
corresponding temperatures are 287 and 350 MeV. The critical temperature 
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size 





statistics 


a [fm] 


T [MeV] 


T/T c 


6 x 20 3 


8.45 


64 + 25 


0.094(1) 


350 


1.17 


6 x 20 3 


8.20 


70 


0.115(1) 


287 


0.96 


16 4 


8.45 


46 


0.094(1) 


113 


0.44 



Table 1: Parameters of our gauge ensembles. 



for the QCD phase transition was determined as T c = 300 MeV for the 
Liischer-Weisz action [T8] . 

We restrict ourselves to configurations with topological charge Q = ±1, 
by considering configurations with a single zero-mode. Thus problems with 
mixing of different zero-modes are avoided. From two larger samples of 
typically 400 configurations we selected 70 configurations with charge Q = 
±1 for j3 = 8.20 and 89 configurations for (3 = 8.45. 

For the deconfined ensemble (f3 = 8.45) the Polyakov loop has a non- 
vanishing expectation value. Due to the Z% symmetry of the gauge ac- 
tion the Polyakov loop can come with three different values for its phase 
(p = 0, ±2-7r/3. The Dirac operator does not have this symmetry and thus 
the eigenvectors and eigenvalues of the Dirac operator will behave differ- 
ently for real, respectively complex Polyakov loop. We found that of the 89 
configurations in the deconfined phase 25 have real Polyakov loop and 64 
have <p = ±27r/3. In the following we will often refer to these two subsets 
as the real and complex sector. 

Finally we also generated a low temperature ensemble on 16 4 lattices at 
[3 = 8.45. It consists of 46 configurations. For all ensembles we use periodic 
boundary conditions for the gauge fields. In Table ^ we summarize the 
parameters of our ensembles. We quote the lattice size, the inverse gauge 
coupling, the statistics, the lattice spacing and the temperature in MeV as 
well as in units of T c . We formally asign also a temperature to the 16 4 
ensemble even though no temporal direction is distinguished by the lattice 
geometry. 

2.2 Observables 

As already outlined in the introduction we use the eigenvectors of the Dirac 
operators to analyze excitations of the QCD vacuum since they provide an 
excellent filter removing the UV fluctuations. Here we use the chirally im- 
proved operator |19| . The chirally improved Dirac operator is a systematic 
expansion of a solution of the Ginsparg- Wilson equation |20| and has very 
good chiral properties. It has been successfully used for quenched spec- 
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troscopy [U with pion masses down to 230 MeV. In [22] it was shown, that 
chirally improved fermions reproduce the analytic result for the zero-mode 
in the background of a discretized instanton very well. In particular for small 
instanton radii it is considerably more accurate than the overlap operator. 

We compute the eigenvectors and eigenvalues of the chirally improved 
Dirac operator using the implicitly restarted Arnoldi method 23 j. For the 
fermions we use periodic boundary conditions for the spatial directions. In 
time direction we implement general boundary conditions that allow for an 
arbitray phase, 

ip(N t + l,x) = exp(i 27rC)^(l,£) • (2-1) 

By Nt we denote the number of lattice points in time direction. The param- 
eter £ can assume values in the interval [0,1], with £ = corresponding to 
periodic temporal boundary conditions and ( = 0.5 giving the anti-periodic 
case. Note that the value £ = 1 gives the same boundary condition as C = 0. 
For all of our configurations as listed in Tabled we computed eigenmodes for 
periodic and anti-periodic temporal boundary conditions. For smaller sub- 
sets from each ensemble we also computed eigenvectors for smaller steps in 
C, such as £ = 0.1, 0.2 ... . More detailed information on these configurations 
is listed in tables below. 

From the eigenvectors ip we construct several gauge invariant observ- 
ables. The simplest observable is the scalar density p(x) of the eigenvectors 
obtained by summing color and Dirac indices c, d for each space-time point 
x separately, 

P(x) = £|V^M 2 . (2.2) 

c.d 

Note that since our eigenvectors ip are normalized one has ^2 X p{x) = 1 
where the sum runs over all lattice points x. From the analytic continuum 
result as well as from lattice studies it is known that the zero-mode traces 
the underlying instanton. Plots of the density p(x) for slices of the lattice 
will be used below to illustrate the behavior of the eigenmodes. We will also 
compare the position of the lumps seen by different eigenmodes ipi,ip2 by 
computing the euclidean 4-distance between their maxima 

d 4 = || x™ ax - x% ax || , (2.3) 

where x™ ax and x™ 0,1 are the lattice points with the largest values of the 
density for tpi, respectively ip2- 
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In order to quantify the localization of the lump seen by the eigenmode 
we use the so-called inverse participation ratio defined as 

/ = F£p(*) 2 , (2.4) 

X 

where V denotes the total number of lattice points. If a mode is localized on 
a single site xq (p(x) = 5 x ,x ) then one finds I = V. For an entirely spread 
out mode (p(x) = 1/V for all x) one has 1 = 1. Thus localized modes have 
large /, while delocalized modes have small values of I. 

Finally we describe an observable to characterize the overlap of lumps 
seen by different eigenmodes. In particular we compute the overlap of the 
support of the two lumps seen by eigenvectors tpi and ip 2 - The support is 
defined as the function 9(x) which has 9{x) = 1 for the lattice points x 
carrying the lump in p and 9{x) = for all other points. An open question 
is how to identify the lump. For a narrow lump the support is small while a 
broad lump has a much larger support. Using the inverse participation ratio 
/ we define a number iV = Int[V/(32/) + 1] and use this to determine the 
support by setting 9{x) = 1 for those N lattice points x which carry the ./V 
largest values of p and 9{x) = for all other points. The normalization of N 
is chosen such that if the lump was a 4-D Gaussian the support would consist 
of the 4-volume carrying the points inside the half-width of the Gaussian. 
Note that N is the size of the support, i.e. N = J2 X 

With our definition of the support of the two lumps (9\ for eigenvector 
ipi, 9 2 for eigenvector ip 2 ) we can define the overlap function C of the two 
lumps as 

C(x) = 9 l {x)9 2 {x) . (2.5) 
From C(x) we compute the relative overlap 

- ivrbil^M- <2 - 6) 

The relative overlap R±2 is a real number ranging from to 1. If the lump in 
pi and the lump in p 2 sit on top of each other (at least the core up to about 
the half-width) and have the same size then R12 = 1. If the two lumps are 
entirely separated one has R12 = 0. Values in between give the ratio of the 
common volume of the two supports to the average size of the two supports. 
Note that even if the two lumps sit on top of each other, but have different 
sizes N±, N2 one has R12 < 1. 
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3 A few words on Kraan - van Baal solutions 



When we present our data in the sections below we will show that the two 
ensembles with finite temperature, i.e. the runs on the 6 x 20 3 lattices, can 
be described by Kraan - van Baal (KvB) solutions. In order to facilitate 
this enterprise let us briefly discuss the structure of the KvB solutions |5| 
and the analytic results for the corresponding zero-mode |13| . 

3.1 Kraan - van Baal solutions 

Kraan van-Baal solutions are classical solutions of the Yang-Mills equations 
on a euclidean cylinder. They are characterized by the Polyakov loop at 
spatial infinity Voo given by 

Voo = lim Pexpf / dtA (t,x)\ , (3.7) 
|2|->oo yJo J 

where j3 is the time extent of the cylinder and the temporal boundary con- 
ditions of the gauge field are periodic, i.e. A^(f3 + t,x) = A^t^x). With a 
suitable constant gauge transformation Voo can be diagonalized and is given 
by (we show the formulae for SU(3)) 

Voo = exp 27rdiag(/ii,^ 2 ,^3)) • (3.8) 

The determinant of Voo has to equal 1 which implies that [X\ + fi2 + ^3 is an 
integer. In order to have a unique labeling of Voo one chooses the constant 
gauge transformation and the phases such that 

^1 + ^2 + ^3 = , 

fJ-i < fJ-2 < M3 < ^4 = 1 + Ml • (3-9) 

In the second line we defined another phase factor \x± which will be useful 
for a compact notation. Using the scale invariance of the classical equations 
the time extent (3 can be set to /3 = 1. Subsequently we use this convention 
and distances are now expressed in units of 0. 

In addition to the phase factors /Uj the KvB solution also depends on 
3 spatial vectors y\ , yi , y% . As already remarked the KvB solution can be 
seen to consist of 3 constituent monopoles and the $ are their positions. As 
for the fi m we define a 4-th vector = y\ for notational convenience. The 
action density of the KvB solutions can be written as 

TrFjL(x) = di^log^(x), 
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Figure 1: Action density of a typical KvB solution for a x,y-slice on 
a logarithmic scale. The parameters of the solution are (^1,^2,^3) = 
(-17, -2, 19)/60, m = (-2, -2, 0), y 2 = (0, 2, 0) and y 3 = (2, -1,0). 



rp(x) = -tr(A 3 A 2 A 1 ) 



cos(27r£) 




Dm Vm+l 

r m +i 



) 




(3.10) 



In these equations r m = \x—y m \, c m = cosh(27Tf m r m ) and s m = sinh(27Tf m r m ) 
Each monopole can be assigned a mass M m , 



This expression for the monopole mass M m implies that the m monopole 
is localized when fi m +i — is large and it is spread out for small /U m _|_i— fi m . 

In order to illustrate the structure of the KvB solutions we show the cor- 
responding action density for an example in Fig.^ The values of the param- 
eters 1 are (/ii, f/Q, ^3) = (— 17, — 2, 19)/60. The positions of the monopoles 
are chosen as y\ = (—2,-2,0), y 2 = (0,2,0) and y 3 = (2,-1,0). We show 
the action density in the x, y-plane on a logarithmic scale. 

3.2 Zero-modes of the Dirac operator 

The eigenvalue problem of the Dirac operator in the background of a KvB 
solution was addressed in ^5] and an explicit expression for the zero-mode 

1 These values are directly taken from the example given in the second reference in |13|. 



M, 



8vr (Mm+l - Mm) ■ 



(3.11) 



S 



was given. In particular the problem was solved for general temporal bound- 
ary conditions (note that we scaled the temporal extent (3 to /3 = 1), 

ip(l + t,x) = exp(i2vrC) tp(t,x) . (3.12) 

The density of the zero-mode is given by 

Pi?) = ElV'^l 2 = -J^ydlf^x) , (3.13) 



with 

f c (x) = —Mx)- 1 (v?\A m - 1 ...AiA a ...A m \w?), (3.14) 




7T 

r.,„ 

and the 2-spinors \v™) and \w™) are given by 

/ sinh(27r(C-^ m )r m )\ , m . / cosh(2vr(C-^„ 
^cosh(27r(C-M m )r m ) > / ' F f ' ^ - sinh(2vr(C-^ 

(3-15) 

The index m which determines the prefactor l/r m , the spinors \v™), \w™) 
and the ordering of the Ai in Eq. (|3.14l) as well as the choice of \i m and 
r m in Eq. (|3.15j) is determined from the interval containing the boundary 
condition parameter, 

fJ-m < C < Mm+1 • (3-16) 

Before we discuss the implications of this selection rule let us denote a 
limiting case of the function f^(x). In the limit of well separated constituent 
monopoles, i.e. large \yi — y%+i\ for all i one finds 

~ sinh(2vr(C-^ m )r m )sinh(27r( / u m+ i-C)r m ) 

td x ) = 2n ^T7o \ • (3-17) 

r m smh(27Tf m r m ) 

This limiting case together with the selection rule Eq. 1)3.16(1 makes explicit 
the following two properties of the zero mode: 

• The zero mode is located on the m th monopole when the boundary 
condition parameter £ is contained in [/j m ,/i m+ i]. The position of the 
lump seen by the zero-mode can jump from y m to y m +i as the value 
of £ crosses ^ m +i- 

• The size of the zero-mode depends in addition to the mass of the 
underlying monopole also on the parameter a, 

a = min{ C - /U m , /U m+ i - C } • (3.18) 

The zero-mode is localized for large values of a and spread-out for 
small a. 
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£=0.1 C=0-5 




C=0.85 



Figure 2: Scalar density of the KvB zero-modes for different values of £ in the 
x, y-plane. The parameters are the same as used in Fig. ^ i.e. (/xi, (12, /U3) = 
(-17, -2, 19)/60, yi = (-2, -2, 0), y 2 = (0, 2, 0) and y 3 = (2, -1, 0). 

In order to demonstrate the features of the zero-modes in Fig. |2 we 
show plots of the scalar density of the zero-modes for the example where 
we discussed the action density in the last subsection and displayed the 
result in Fig. ^ Our example has fj,2, A*3, Ha) = (— 17, — 2, 19, 43)/60. 
This implies that for £ £ [—2/60, 19/60] the zero mode is located on y%, for 
C G [19/60,43/60] it is located on y 3 and for C G [43/60,58/60] it is located 
on y\. Note that we made use of the equivalence £ = Q — 1 in order to treat 
the interval [^1,^2] = [-17/60,-2/60]. 

The plots show nicely that indeed the position of the lump seen by the 
zero-mode can change its position. In the following sections we now make 
use of this property and use the zero-modes of the lattice Dirac operator 
with different boundary conditions as a device for detecting different lumps 
that together build up a topological excitation of charge Q = ±1. 
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4 Finite Temperature - results above T c 



4.1 Overview over the configurations 

We remarked in Section 2 that the ensemble in the deconfined phase (6 x 20 3 , 
(3 = 8.45) consists of 89 configurations, 25 with real Polyakov loop and 
64 with complex Polyakov loop. For a subset of 10 configurations we did 
not only compute the zero-mode with periodic and anti-periodic tempo- 
ral boundary conditions but for a general phase exp(i27r£) at the temporal 
boundary. In particular we studied Q = 0.0, 0.1, 0.2, 0.3, 1/3, 0.4, 0.5, 0.6, 
2/3, 0.7, 0.8, 0.9. The cases £ = 0.0 and Q = 0.5 correspond to periodic 
respectively anti-periodic boundary conditions. For this subset Table |2] lists 
some basic properties. In particular we list the configuration number, the 
phase of the Polyakov loop, the inverse participation ratio for the periodic 
and the anti-periodic zero-mode, the position of the maxima for periodic 
and anti-periodic b.c, and finally the 4-distance between the two maxima. 

A few basic properties of the configurations above T c can be immedi- 
ately read off from Table [21 For all configurations with real Polyakov loop 
{ip = 0) the inverse participation ratio of the zero-mode increases when 



conf. 


f 


Ip 


Ia 


(t 


x,y 


Z)p 


(t,x,y 


„\max 

z )a 


d±/a 


6 





2.19 


43.16 


5 


10 


9 


8 


3 


4 


17 


9 


10.24 


7 


-2tt/3 


22.81 


4.71 


3 


3 


3 


11 


1 


4 


2 


9 


3.16 


9 


+2vr/3 


23.70 


5.93 


5 


8 


18 


17 


5 


8 


18 


17 


0.00 


10 


+2vr/3 


40.27 


10.04 


1 


11 


13 


14 


1 


11 


13 


14 


0.00 


19 





2.29 


20.20 


5 


11 


9 


8 


4 


20 


3 


8 


10.86 


38 





2.38 


31.00 


5 


12 


9 


8 


5 


12 
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8 


0.00 


49 


+2vr/3 


34.03 


15.67 


6 


11 


19 


2 


6 


11 


19 


2 


0.00 


87 





1.90 


47.76 


6 


18 


13 


9 


6 


18 


13 


9 


0.00 


101 





2.87 


41.69 


5 


12 


18 


10 


1 


19 


3 


1 


12.60 


383 





2.42 


42.86 


6 


19 


13 


11 


2 


3 


11 


4 


8.54 



Table 2: Parameters of the subensemble in the deconfined phase (6 x 20 3 , 
(3 = 8.45) which we studied with finely spaced values for the boundary 
parameter £. For these 10 configurations we quote the configuration number, 
the phase of the Polyakov loop, the inverse participation ratios Ip and I a for 
the zero-modes with periodic, respectively anti-periodic b.c, the position of 
the maxima for the corresponding zero-modes and the 4-distance d± between 
these two maxima in units of the lattice spacing. 
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switching from periodic to anti-periodic boundary conditions, i.e. the mode 
becomes more localized. For the configurations with complex Polyakov loop 
{ip = ±27r/3) the situation is reversed and the state becomes more delocal- 
ized when switching from periodic to anti-periodic b.c. Furthermore only 
configurations with real Polyakov loop allow for a large distance between the 
lumps seen by the periodic and anti-periodic zero mode. In the following 
subsections we will elaborate further on these observations. 

4.2 Example for the generic behavior 

In this subsection we show series of plots of the scalar density for a config- 
uration where we computed with the zero-modes with finely spaced values 
of Q. Before we do so let us discuss what can be expected from the analytic 
KvB result. 

In the deconfined phase the Polyakov loop can have the three phases 
if = 0, ±27r/3. Let us first discuss the case of real Polyakov loop. For this 
case one finds for the phase factors \ii entering the KvB solutions, 

(/ii,M2,/i 3 ,/i4) = (0,0,0,1) . (4.19) 

Applying the rule that the zero mode with boundary factor £ is located at 
position y m when ( £ [fi m , fi m+ i] we find: For all £ ^ the zero mode 
is located on y^. The only exception is £ = where the zero mode is 
located on y\ and yi ■ Note that £ = is equivalent to ( = 1 such that also 
for ( = 1 this second possibility holds. The mode at ( = 0, 1 should be 
very much delocalized, while the other mode seen for £ ^ is expected to 
show a changing degree of localization. In particular one expects a rather 
delocalized mode near the endpoints of the interval £ = 0, 1 while it is 
more localized in the middle of the interval, reaching a maximum of the 
localization at £ = 0.5. 

The plots we will show below are for our configuration 383 from the 
6 x 20 3 ,/3 = 8.45 ensemble. In Table El we show the inverse participation 
ratio and the position of the maximum of the scalar density for all values 
of C- From the table it is obvious that the mode is located at the same 
position for all ( ^ 0, 1 and at a different position for £ = 0, 1. The distance 
between the two lumps is d±/a = 8.54 in lattice units or in fermi d± = 0.80 
fm. The inverse participation ratio is largest near £ = 0.5, i.e. the mode is 
most localized for anti-periodic boundary conditions. At £ = 0, 1 the inverse 
participation ratio is small and the mode is very much spread out. 

Let us now look at the plots of the scalar density for configuration 383. 
In Fig. Elwe show the scalar density on x, y-slices of the lattice. In particular 
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0.0 


2.42 


6 


9 


13 


11 
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45.32 
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9.15 
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13 


11 
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0.8 


17.00 
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13 


11 4 


0.3 


19.05 
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13 


11 
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0.9 


3.40 


2 


13 


11 4 


0.4 


31.80 


2 


13 


11 


4 


1.0 


2.42 


6 


9 


13 11 


0.5 


42.86 


2 


13 


11 


4 













Table 3: Inverse participation ratio and position of the maximum of the 
scalar density as a function of the boundary condition parameter £. The 
data are for configuration 383 of the deconfined ensemble (6 x 20 3 , f3 = 8.45). 
Plots of the corresponding scalar density are shown in Fig. El 

we show the slice at t, z = 6, 11 in the left column of plots and the slice at 
t, z = 2, 4 in the right column of plots. Thus in the left column we show 
the slice where the scalar density at £ = 0, 1 has its peak while in the right 
column the slice for the maximum at all other values of £ is shown. The 
reason for showing two different slices is that typically for a thermalized 
configuration the positions of the two maxima do not fall in a common 
coordinate plane. In the plots we use values of £ = 0, 0.2 and 0.5. Note that 
the scale for the right column differs by a factor of 10 from the scale in the 
left column. 

For C = there is a clearly visible lump in the 6,11-slice (left column, top 
plot of Fig. EJ) , while the 2,4-slice is essentially flat (right column, top plot). 
As we increase £ to £ = 0.2 the peak in the 6,11-slice has vanished, but a 
small lump in the 2,4-slice has emerged (second row in Fig. OJ). This latter 
lump starts to grow as we increase £ and reaches its maximum in the last 
row of Fig. 01 at £ = 0.5. As we increase ( further the lump starts to shrink 
again, and essentially runs through the series of plots in reverse order. The 
series of plots shown in Fig. 01 nicely demonstrates that the zero-mode for 
configuration 383 behaves exactly like the zero-mode for a KvB solution in 
a deconfined configuration with real Polyakov loop. 

Let us now look at the other two possibilities for the phase of the 
Polyakov loop. Phases of ip = 2ir/3 and ip = —2tt/3 correspond to the 
following two sets of jii, 

(/ix,^,^,^) = (-2/3,1/3,1/3,1/3), (4.20) 
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Figure 3: Slices of the scalar density for 6 x 20 3 , P = 8.45 and real Polyakov 
loop (configuration 383). We show x,y-slices at t = 6, z = 11 (left column) 
and at t = 2, z = 4 (right column). The values for £ are ( = 0, 0.2, 0.5 (from 
top to bottom). Note the different scale for the l.h.s. and r.h.s. plots. 



for ip = 27r/3 and 

( W , M2)M3) M4) = (-1/3,-1/3,2/3,2/3), (4.21) 

for tp = — 27r/3. Thus for tp = 2n/3 the zero-mode has to be located at the 
same position for all £ ^ 1/3 and only for £ G [1/3, 1/3] (in other words for 
£ = 1/3) the lump can sit at a different position. For (p = 2tt/3 the zero- 
mode has to be located at the same position for all £ ^ 2/3, but can have 
a different location for £ = 2/3. This discussion shows that the situation 
for (p = ±27r/3 is equivalent to the <p = case and only the interval [0, 1] is 
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shifted, such that the position of a possible jump of the location is shifted 
to C = 1/3, respectively £ = 2/3. 

We find that all our configurations in the deconfined ensemble obey 
the general pattern we have discussed. Certainly it is also clear that our 
thermalized configurations are not classical solutions and quantum effects 
play a role. For example we find for the configuration shown in Fig. [31 that 
the peak in the 6, 11-slice is already visible at £ = 0.9. For an unperturbed 
KvB solution it should be visible only at exactly ( = 1 (equivalent to £ = 0). 
The quantum effects seem to slightly twist the Polyakov loop such that it 
is not exactly an element of the center. This makes the critical values of 
£ = 0, 1/3,2/3 where the lump can change its position less strict than for 
the continuum solution. 

4.3 Results from the whole sample 

In the last subsection we have discussed an example for the generic behav- 
ior of the KvB zero-modes in the deconfined phase and remarked that all 
configurations we looked at show the same general pattern. Here we now 
present results for the whole set of configurations listed in Table [21 where 
we used finely spaced C and also for the even larger set of 89 configurations 
where we only compared periodic and anti-periodic boundary conditions. 

Let us begin with looking at the behavior of the inverse participation 
ratio, i.e. the localization property. The discussion in the last section and 
the example shown there suggest that for all configurations in the real sector 
the inverse participation ratio has a maximum near £ = 0.5 (the zero- mode is 
most localized there) and a minimum near £ = 0, 1. For configurations with 
complex Polyakov loop the whole curve is shifted such that the minimum 
is at £ = 1/3 for ip = 2tt/3 and at £ = 2/3 for <p = —2tt/3. Note that this 
behavior is independent from a possible change of the position of the peak. 

In Fig. 0] we show for all configurations in the deconfined phase where 
we did runs with small steps for £ how the localization changes as a function 
of £. The l.h.s. plot gives the results for configurations with real Polyakov 
loop, while the r.h.s. plot is for complex Polyakov loop. 

The plots show that indeed for all these configurations the predicted 
pattern holds. The configurations with real Polyakov loop (l.h.s. plot) have 
a maximum of I near £ = 0.5 and minima near and 1. For configurations 
with complex Polyakov loop (r.h.s. plot) we find the curve shifted such that 
the minimum is at £ = 1/3 for phase </? = 27r/3 (configurations 9, 10, 49 in 
the r.h.s. plot) and at C = 2/3 for phase (p = —2ir/3 (configuration 7). 

From the above discussions of the properties of KvB zero-modes and 
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Figure 4: Inverse participation ratio / as a function of (. We show the 
results for the deconfined phase (6 x 20 3 , (3 = 8.45), displaying configurations 
with real Polyakov loop in the l.h.s. plot and configurations with complex 
Polyakov loop in the r.h.s. plot. 

the examples we gave one can derive another prediction. This prediction 
links the distance between the positions of the lumps in the zero-mode with 
periodic and anti-periodic boundary conditions to the phase of the Polyakov 
loop: For real Polyakov loop the periodic mode with £ = can sit in either 
the interval [0,0] or in the interval [0, 1], while the periodic zero-mode with 
C = 0.5 is only contained in [0, 1]. Thus for real Polyakov loop the periodic 
and the anti-periodic zero-mode can be located at different positions. The 
situation is different for complex Polyakov loop. E.g. for ip = —2n/3 both 
C = and C = 0.5 are contained in the interval [—1/3, 2/3] and the periodic 
as well as the anti-periodic zero-mode are expected to be located at the same 
position. 

In Fig. |S]we test this property and in the l.h.s. plot show histograms of 
the distance between the peaks in the scalar density when comparing pe- 
riodic to anti-periodic temporal boundary conditions. From the histogram 
it is obvious that indeed all configurations with large distances between 
the peaks, i.e. distances larger than 5 lattice spacings (0.47 fm) have real 
Polyakov loop. About a third of the configurations with complex Polyakov 
loop also show small but non-zero distances (up to distance 5 for a single 
configuration) between the periodic and the anti-periodic lump. We believe 
that this is due to quantum fluctuations superimposed on the lumps. Note 
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Figure 5: Histograms for the distance Rap between the peaks in the scalar 
density when comparing periodic to anti-periodic b.c. (l.h.s. plot) and for 
the corresponding relative overlap Rap (r.h.s. plot). The data are for the 
deconfined ensemble on 6 x 20 3 lattices at f3 = 8.45. 



that also for real Polyakov loop the KvB zero modes can sit at the same 
position a case which is realized by about half of the zero modes in the real 
sector. The plot on the r.h.s. of Fig. \5\ shows histograms for the relative 
overlap Rap of the periodic and the anti-periodic zero mode. These results 
corroborate the findings of the histograms for the distance since the config- 
urations which have no or only small overlap are all in the real sector. Note 
that a relative overlap smaller than 1 can also come from two modes which 
sit at the same position but have different localization. This explains why 
also configurations in the complex sector have relative overlap smaller than 
1 (compare also Fig. Eland the corresponding discussion below). 

Finally in Fig. El we show a scatter plot of the values of the inverse par- 
ticipation ratios Ip,Ia for periodic and anti-periodic boundary conditions. 
For the real sector the mode is delocalized and Ip is small. Thus one ex- 
pects that in the Ip versus 1a scatter plot the data fall on a horizontal line. 
Our numerical results clearly show this behavior. For the complex sector 
the periodic mode has a scale parameter a ~ 1/3 while the anti-periodic 
zero mode has a ~ 1/6. Thus one expects Ip > I a for the complex sector 
and again our data nicely confirm this prediction. The plot also contains 
numbers of Ip and I a that were generated from the analytic solution for 
complex Polyakov loop. We varied the positions yi of the monopoles thus 
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Figure 6: Scatter plot for the inverse participation ratio of the zero mode 
with periodic and anti-periodic boundary conditions. The data are for the 
deconfined ensemble on 6 x 20 3 lattices at (3 = 8.45. We use filled circles for 
configurations with real Polyakov loop and open triangles for the complex 
sector. We also include numbers (crosses connected with straight lines to 
guide the eye) that were generated from the analytic solution. 

obtaining different values for the inverse participation ratio. These data are 
shown as crosses in the plot and we connect them to guide the eye. Our 
numerical results in the complex sector follow this analytic curve quite well. 
Only for very strongly localized modes {Ip > 100) we see a slight deviation. 
We remark that on the lattice such strongly localized states are resolved by 
only a few lattice points such that cut-off effects become important. 

4.4 Results for the spectral gap 

The chiral condensate is related to the spectral density /9 spe c(A) of the Dirac 
operator near the origin through the Banks-Casher formula [5], 

(^V) = - Jim ^Ps P ec(0). (4.22) 

For the confining phase where chiral symmetry is broken the spectral den- 
sity extends all the way to the origin. In the deconfined phase where chiral 
symmetry is restored a gap opens up in the spectral density near the ori- 
gin. /O sp ec(0) vanishes and so does the condensate. The opening up of the 



18 



Configurations with real Polyakov loop Configurations with complex Polyakov loop 




Figure 7: Size of the spectral gap as a function of £. We show the re- 
sults for the deconfined phase, displaying configurations with real Polyakov 
loop in the l.h.s. plot and configurations with complex Polyakov loop in the 
r.h.s. plot. 

spectral gap is believed to come from an arrangement of instantons and 
anti-instantons to tightly bound molecules [5]. 

An interesting question is whether the chiral condensate vanishes at the 
same temperature for all sectors of the Polyakov loop. Results [21] for 
staggered fermions with anti-periodic temporal b.c. seemed to indicate that 
the chiral condensate vanishes at different temperatures for real and complex 
Polyakov loop, and attempts to understand this phenomenon can be found in 
|25j . A more recent study |18j analyzing directly the spectral gap (also with 
anti-periodic temporal b.c.) came to a different conclusion: The spectral 
gap appears at the same critical temperature for all sectors of the Polyakov 
loop and this temperature coincides with the deconfinement transition. Here 
we further corroborate the scenario of a single transition by discussing the 
spectral gap for different boundary conditions for the Dirac operator. 

In Fig. [7| we show the spectral gap as a function of the boundary con- 
dition parameter £ for the 10 configurations in the deconfined phase where 
spectra with finely spaced £ were computed. The l.h.s. plot shows the 6 
configurations with real Polyakov loop, while the r.h.s plot gives data for 
the configurations with complex Polyakov loop. The spectral gap is defined 
as the size of the imaginary part of the smallest complex eigenvalue. 

The plots show that as for the inverse participation ratio the results for 
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the gap from configurations with <p = ±2n/3 for the Polyakov loop can be 
obtained from the result for <p = by shifting the £ axis by amounts of 
1/3 and 2/3 respectively. The results in the three sectors are completely 
equivalent (with only the data for configuration 38 showing a somewhat 
irregular behavior) and can be transformed into each other through different 
boundary conditions for the fermions. This symmetry shows that there is 
no reason to expect a different critical temperature for chiral symmetry 
restoration in different sectors of the Polyakov loop. 

It is interesting to note that in the real sector for £ = and in the com- 
plex sectors for £ = 1/3, respectively £ = 2/3 the spectral gap becomes very 
small. This could indicate that for these boundary conditions chiral symme- 
try remains broken also in the deconfined phase. It would be interesting to 
analyze the chiral condensate in the deconfined phase using these boundary 
conditions for the fermions. 

5 Finite Temperature - results below T c 

5.1 Overview over the configurations 

In the last section we have studied the behavior of the configurations in the 
deconfined phase. We have found that qualitatively our results match the 
characteristics of the KvB solutions very well. In this section we now turn 
to the ensemble below T c . Below T c the Polyakov loop vanishes which leaves 
only the unique possibility, 

(Ati,M2,M3,M4) = (-1/3,0,1/3,2/3), (5.23) 

for the phase factors /ij. One expects that the zero-modes can change their 
position when £ crosses the values of 1/3 and 2/3. Thus below T c the zero- 
mode can change its position more often and for these configurations we 
now use a finer spacing of 

From the complete ensemble below T c (compare Table ^) we computed 
for 10 configurations eigenvectors at all values of ( between £ = and 
£ = 1 in steps of 0.05. In Table |1] we list some basic properties of these 
configurations. In particular we list the configuration number, the inverse 
participation ratios Ip and I a for the zero-modes with periodic, respectively 
anti-periodic b.c, the position of the maxima for the corresponding zero- 
modes and the 4-distance between these two maxima in lattice units. 
Finally we list the number of different lumps visited by the mode in a com- 
plete cycle through £. We count a lump as independent lump when its 



20 



COI1I. 


Tr, 


T A 




t, x, y 


Z) p 




(t 


x -> Vi 


z )a 


ax 




~fr hitrinc 

■j-h ILllIlUo 


13 


4.67 


13.41 


5 


11 


16 


11 


2 


1 


9 


6 


13.52 


4 


18 


7.32 


10.73 


2 


6 


20 


3 


6 


16 


15 


12 


14.49 


3 


92 


38.52 


5.73 


1 


11 


13 


13 


4 


8 


16 


11 


5.56 


2 


109 


5.56 


15.93 


1 


14 


18 


18 


2 


18 


5 


3 


9.53 


4 


123 


4.59 


32.21 


5 


18 


19 


11 


4 


13 


18 


1 


11.26 


2 


125 


4.13 


14.72 


5 


1 


16 


9 


5 


14 


16 


18 


11.40 


3 


153 


22.98 


10.32 


5 


17 


18 


19 


6 


18 


17 


2 


3.46 


2 


210 


2.48 


28.96 


1 


1 


8 


4 


2 


2 


1 


13 


11.48 


3 


215 


7.23 


3.63 


1 


15 


20 


14 


5 


11 


17 


6 


9.64 


3 


266 


8.15 


40.01 


1 


9 


19 


20 


2 


9 


6 


6 


9.27 


3 



Table 4: Parameters of the subensemble below T c (6 x 20 3 ,/3 = 8.20) which 
we studied with finely spaced values for the boundary parameter £. For these 
10 configurations we list the configuration number, the inverse participation 
ratios Ip and I a for the zero-modes with periodic, respectively anti-periodic 
b.c, the position of the maxima for the corresponding zero-modes, the 4- 
distance between these two maxima and the number of different lumps 
visited by the mode in a complete cycle through (. 

maximum remains at the same position for at least 3 subsequent values 
of C, i.e. dominates the zero mode at least for a ("-interval of length 0.1. 
Note that these lumps are visited by the zero-mode one after the other and 
typically one of them dominates the zero-mode. 

A first interesting finding can already be read off from the numbers in 
the table. As discussed, below T c the zero mode has three possible values of 
( where it can change position, ( = 0, 1/3, 2/3. Thus the zero-mode can visit 
up to 3 different positions. Also a number less than 3 is possible since the 
coordinate vectors yi of the monopoles can coincide. Most of the eigenmodes 
do indeed see 2 or 3 lumps but for two of the configurations (13 and 109) we 
found 4 lumps. This observation indicates that for our ensemble below T c 
the zero-modes do not follow the predictions for KvB zero-modes as closely 
as their counterparts above T c . Reasons for that may be quantum effects 
which drive the trace of the Polyakov loop away from zero. We also remark 
that the lattice below T c is coarser than the lattice above T c . 
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5.2 Example for the generic behavior 

As for the ensemble in the deconfined phase we begin the presentation of 
our results with showing slices of the scalar density for a generic example. 
As discussed in the last subsection, below T c the zero-mode is expected to 
visit up to three different lumps in a complete cycle through the boundary 
condition. Configuration 125 is an example for such a behavior. In Table 
we list some basic properties of this configuration for all values of £ we 
studied. 

Table |S] shows that the zero- mode of configuration 125 visits three differ- 
ent lumps located at (t, x, y, z) = (5,13,6,9), (2,6,13,19) and (5,6,6,18). These 
three lumps are well separated from each other with mutual 4-distances of 
e?4 = 14.38, d± = 11.40 and d^ = 7.68 in units of the lattice spacing. In 
fermi this is d^ = 1.65 fm, 1.31 fm and 0.88 fm. The most prominent lump 
is the one at (i, x, y, z) = (5, 6, 6, 18) which is seen by the zero mode at all £ 
between £ = 0.4 and £ = 0.85. We remark that this lump is also the most lo- 
calized one, i.e. reaches the largest values for the inverse participation ratio. 
This is as expected for the KvB zero-modes, where the mass of the under- 
lying monopole is proportional to /i m +i — [i m (compare Ea. l3.llj) . Thus the 
lump which occupies the largest interval in £ is also the most localized one. 
This pattern also holds for all other configurations listed in Table 0J It is 
interesting to note that the lump at (t, x, y, z) = (5, 6, 6, 18) already briefly 
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Table 5: Inverse participation ratio and position of the maximum of the 
scalar density as a function of the boundary condition parameter £. The 
data are for configuration 125 of the confined ensemble (6 x 20 3 , (5 = 8.20). 
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Figure 8: Slices of the scalar density for 6 x 20 3 ,/? = 8.20, configuration 
125. We show x,y-slices at t = 5, z = 9 (left column), at t = 2, z = 19 
(center column) and t = 5, z = 18 (right column). The values for £ are 
C = 0.05,0.3,0.65 (from top to bottom). 

appears at £ = 0.1,0.15 before the lump at (t,x,y,z) = (2,6,13,19) takes 
over. The reason for that is that the lump at (2,6,13,19) starts to grow only 
slowly such that for the small window in £ the remains of the tallest lump 
at (5,6,6,18) dominate. 

Let us now look at plots for the scalar density for this configuration. In 
Fig. [H] we show x,y-slices for the scalar density taken through the maxima 
of the three lumps, i.e. at t = 5, z = 9 (left column), at t = 2, z = 19 
(center column) and t = 5, z = 18 (right column). The values for £ are 
C = 0.05, 0.3, 0.65. These are the values where the dominating lumps have 
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Figure 9: Slices of the scalar density for 6 x 20 3 , (3 = 8.20, configuration 125. 
We show t, x-slices at y = 6, z = 9 for £ = 0.05 (left plot), at y = 13, z = 19 
for Q = 0.3 (center plot) and y = 6, z = 18 for £ = 0.65 (right plot). These 
plots correspond to the three dominant lumps also shown in Figs. |H1 



their maximum. In between these extremal values one has intermediate 
situations with a receding and an advancing lump. Note that the plots in 
the r.h.s. column have a different scale and the corresponding lump is about 
twice as tall as the other two lumps. 

It is interesting to inspect also other slices, in particular a slice containing 
the time direction. For an unperturbed KvB solution one expects that 
the lump extends over all of the time direction and that the height of the 
lump displays an oscillating behavior in time. For the action density of 
cooled SU(2) lattice configurations with twisted boundary conditions such 
a behavior was observed in |llj . 

In Fig. |2lwe show t, x-slices of the scalar density again for configuration 
125 which was already used for the plots in Fig. |HJ In particular we show the 
slices through the 3 lumps at their dominant value of £. They were taken 
at y = 6, z = 9 for £ = 0.05 (left column), at y = 13, z = 19 for £ = 0.3 
(center column) and y = 6, z = 18 for £ = 0.65 (right column). The plots 
show that all three lumps are stretched along the i-axis. A similar behavior 
was also observed for other configurations, both below and above T c . 

5.3 Results from the whole sample 

As for the deconfined phase, also for the ensemble below T c we now study 
observables evaluated for the whole ensemble in order to establish properties 
of the zero-modes beyond a demonstration in a single example. 

We start with plots for the inverse participation ratio as a function of (. 
In Fig. E3 w e show such figures for the 10 configurations with finely spaced 
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Figure 10: Inverse participation ratio / as a function of £ for the 6 x 20 3 , 
P = 8.20 ensemble. In order to avoid overcrowded figures we show the data 
for configurations 13, 18, 92, 109 and 123 on the l.h.s., and configurations 
125, 153, 210, 215, 266 on the r.h.s. 

C listed in Table |1J To avoid overcrowding of the plots we split the figure 
into two plots and show the data for configurations 13, 18, 92, 109 and 123 
on the l.h.s., and configurations 125, 153, 210, 215, 266 on the r.h.s. 

From the plots it is obvious that the localization of the zero-mode fluc- 
tuates strongly during a complete cycle through the boundary condition. 
Typically there is a single pronounced maximum in the inverse participa- 
tion ratio and one or two local maxima. The tallest peak occupies also the 
largest interval of £ in accordance with the mass formula Eq. (|3.11j) . On the 
other hand it is obvious from the plots that the minima are not necessarily 
located at exactly £ = 0, £ = 1/3 and C = 2/3 as would be the case for an 
unperturbed KvB zero-mode. However, one should not forget that our con- 
figurations are taken from a thermalized ensemble containing all quantum 
effects. It is not yet known [26^ in which way quantum effects modify the 
classical KvB solutions. 

In Fig.^Jwe show histograms for the distance between the lump in the 
periodic and anti-periodic zero-modes in the l.h.s. plot. On the r.h.s. we show 
a histogram for the relative overlap Rap between these two lumps. These 
two histograms contain the data for all 70 configurations in the ensemble. 

The histogram on the l.h.s. shows that for half of the configurations the 
zero-mode for periodic boundary conditions and its anti-periodic counterpart 
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Figure 11: Histograms for the distance d± between the peaks in the scalar 
density when comparing periodic to anti-periodic b.c. (l.h.s. plot) and for 
the relative overlap Rap between the two zero- modes (r.h.s. plot). The data 
are for the ensemble in the confining phase on 6 x 20 3 lattices at /3 = 8.20. 

sit on lumps which are more than 5 lattice spacings apart, i.e. at least 0.69 
fm. The plot on the r.h.s. supports this finding by illustrating that typically 
the two lumps do not overlap substantially. Note that even if the lumps 
are located at the same position their relative overlap can be smaller than 
1 when they differ in size. Both plots in Fig. ^2 show that configurations 
where the zero-mode changes its location and size are not singular events 
but make up a large portion of the ensemble. 

6 Results on the torus 

6.1 Overview over the configurations 

Having successfully identified the substructure predicted for the KvB zero- 
modes in ensembles with high temperature it is natural to apply the same 
techniques to configurations at low temperature. In this section we perform 
such an exploratory study and analyze our ensemble on 16 4 lattices with 
(3 = 8.45. The difference to the high temperature configurations is that 
the lattice geometry does not single out a time direction. We break the 
euclidean symmetry explicitly and use our generalized fermionic boundary 
condition for the 1-direction, while for the other directions we implement 
periodic boundary conditions. We remark that on the torus no analytic 
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results are available to compare with. 

As for the cases with temperature we compute eigenvectors of the Dirac 
operator with different values of the boundary condition parameter £. For 
all 46 configurations of the ensemble (compare Table^) we use £ = and ( = 
0.5, i.e. periodic and anti-periodic boundary conditions. For a subensemble 
of 10 configurations we computed eigenvectors for all values of C between 
and 1 in steps of 0.1. 

For these 10 configurations Table El shows some basic properties. In 
particular we list the configuration number, the inverse participation ratio 
for the periodic and anti-periodic zero-modes, the positions of the maxima 
in the corresponding scalar densities and the difference between these two 
peaks. In addition we display the number of lumps visited by the zero-mode 
in a complete cycle through the phase at the boundary. Here we count a 
lump as independent if it appears at two subsequent values of (, i.e. is visible 
in a ^-interval of at least length 0.1. This is the same criterion used for the 
6 x 20 3 ,/3 = 8.20 ensemble. 

The table makes obvious that also for the ensemble on the torus we do 
find configurations where the zero-mode is located at different positions for 
different values of The 4-distances between the periodic and the anti- 
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Table 6: Parameters of our ensemble on the torus (16 , j3 = 8.45). We list 
the configuration number, the inverse participation ratios Ip and I a for the 
zero- modes with periodic, respectively anti-periodic b.c, the position of the 
maxima for the corresponding zero-modes, the 4-distance d^ between these 
two maxima and the number of different lumps visited by the mode in a 
complete cycle of (. 
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periodic lumps reach values up to 11 in lattice units which is up to 1.03 
fm. When discussing the complete ensemble we will find that about half 
of the configurations have a sizeable distance between the two peaks. The 
zero-modes visit between 1 and 3 different lumps in a complete cycle of 
£• We remark that although for configuration 561 we have ^ we list 
only a single lump in the last column. The reason is that only for £ = 
the zero-mode is located at a different position and thus the criterion of 
two subsequent appearances is not fulfilled. The pattern observed for this 
configuration is similar to the pattern of a configuration in the deconfined 
phase with real Polyakov loop. 

6.2 Example for the generic behavior 

Let us now look at figures of the scalar density for a typical configuration. 
In particular we will show plots for configuration 32 of our ensemble. Before 
presenting these figures we discuss Table where we list some properties of 
this configuration as a function of £. In the left half of the table we show £, 
the inverse participation ratio and the position of the maximum. The zero- 
mode changes its position at £ = 0.3 and £ = 0.9. The localization shows 
some mild changes with reaching local maxima at £ = 0.5 and £ = 0.1, 
i.e. approximately in the middle of the ^-intervals for the two lumps. 

In the right half of the table we show the data from a small consistency 
check which we performed for 4 of the 10 configurations. We switched the 
direction where we apply the boundary condition from the 1-direction to the 
2-direction. Since our lattice does not single out a particular direction any 
direction for applying the boundary condition serves equally well and the 
results should be consistent. The data in the right half of the table show that 
this is indeed the case. Also with the non-trivial boundary condition applied 
in the 2-direction we see the same two lumps at (5,9,5,8) and (15,5,8,15). 
They appear, however, at different values of ( which is not surprising since 
a single gauge configuration does not respect rotational symmetry. Also 
for the other 3 configurations where we experimented with changing the 
direction we single out by the boundary condition, we find that essentially 
the same lumps are seen, again at slightly different values of £. 

In Fig. ^1 we show £2, X3-slices Q f the scalar density for configuration 32 
(in the plots we again label these two axes as x and y). The slices were taken 
through the maxima of the two peaks, i.e. at xi,X4 = 5,8 (left column of 
plots), respectively at xi,x$ = 15,15 (right column of plots). We use values 
of £ = 0, 0.3, 0.5 showing the two extremal situations where only one of the 
lumps dominates (£ = 0, top) and ( = 0.5, bottom) and an intermediate 
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Table 7: Inverse participation ratio and position of the maximum of the 
scalar density as a function of the boundary condition parameter £. The 
data are for configuration 32 of the ensemble on the torus (16 4 ,/3 = 8.45). 
The data in the left half of the table were obtained with using the nontrivial 
boundary condition in 1-direction, while the data in the right half is for 
2-direction. The corresponding plots of the scalar density are shown in Fig. 

□21 

situation (( = 0.3, middle). Note the different scale for the two columns 
of plots. The data we show are from the run with the boundary phase 
attached in 1-direction. The plots show clearly that the zero-mode changes 
its location with The same pattern with certain values of £ where a 
single lump dominates and intermediate values with an advancing and a 
retreating lump is also seen for the other configurations which visit two 
or three lumps. Equivalent plots taken for the zero-modes computed with 
boundary conditions in 2-direction show similar behavior. The lumps sit 
at the same positions and even have similar shape, such as e.g. the slightly 
elongated form along the x-axis in the plots in the left columns. 

6.3 Results from the whole sample 

Let us now try to analyze the common features of all configurations in the 
ensemble. We start with discussing the inverse participation ratio as a func- 
tion of £. In Fig.EHwe show two such figures where we again distribute the 
10 configurations among two plots, with the l.h.s. plot containing configu- 
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Figure 12: Slices of the scalar density for 16 4 ,/3 = 8.45, configuration 32. 
We show X2, 23-slices (in the plots we again label these two axes as x and 
y) at x\ = 5, X4 = 8 (left column) and at x\ = 15, 24 = 15 (right column). 
The values for Q are (from top to bottom) Q = 0,0.3,0.5. Note the different 
scale for the l.h.s. and r.h.s. plots. 



rations 32, 64, 82 and 93 and configurations 103, 351, 354, 522, 561 on the 
r.h.s. We remark that configuration 61 which is also listed in Table H3 has 
values of / above 50 for all £ such that this defect-like lump does not show 
up for the plot-range we have chosen. 

The plots show that some of the configurations display quite a sizeable 
change in their localization as a function of Q. Four of the zero-modes (64, 
93, 103, 351) show only relatively small changes of their inverse partici- 
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Figure 13: Inverse participation ratio I as a function of £ for the 16 4 , = 
8.45 ensemble. In order to avoid overcrowded figures we show the data for 
configurations 32, 64, 82 and 93 on the l.h.s., and configurations 103, 351, 
354, 522, 561 on the r.h.s. 

pation ratio. An inspection of Table shows that this behavior seems to 
be correlated with neither e?4 nor with the number of lumps visited by the 
zero-mode. 

Let us now address the question whether changing its location is a com- 
mon property of the zero-modes for our ensemble on the torus. We do this 
by again analyzing histograms for the distance d^ between the maxima of 
the periodic and anti-periodic scalar density and histograms for the overlap 
between the two lumps. In Fig. El we show such plots with the histogram 
for the distance on the l.h.s. and the histogram for the overlap on the r.h.s. 
The histograms were taken from all 46 configurations in our ensemble on 
the torus. 

The histogram for the distance shows that indeed configurations are 
quite abundant where the zero-mode changes its position by a sizeable 
amount, larger than a simple deformation by a fluctuation. In particular 
about half of the configurations show a d& of at least 6 in lattice units, 
i.e. the centers of the lumps are at least 0.56 fm apart. The distance d^ 
reaches values up to 13 in lattice units, i.e. 1.22 fm. The histogram for 
the relative overlap shows that even most of the zero-modes that do not 
change their location suffer some deformation when changing the boundary 
condition, such that the relative overlap is pushed below 1. 
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Figure 14: Histograms for the distance between the peaks in the scalar 
density when comparing periodic to anti-periodic b.c. (l.h.s. plot) and for 
the relative overlap (r.h.s. plot). The data are for the ensemble on the torus 
(16 4 , /? = 8.45). 




Figure 15: Scatter plot for the inverse participation ratio of the zero mode 
with periodic and anti-periodic boundary condition (Ip respectively I a)- 
The data (triangles) are for the ensemble on the torus (16 4 , = 8.45). We 
display the line with Ip = Ia dashed line. 
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Finally in Fig. ^] we show a scatter plot of the periodic inverse partic- 
ipation ratio Ip versus its anti-periodic counterpart I a- For the deconfined 
phase the equivalent plot shown in Fig. Hollas revealed an asymmetry between 
periodic and anti-periodic boundary conditions which is closely related to 
the critical values of £ where the zero-mode can change its location. For 
the ensemble on the torus our Fig. [H)\ does not reveal such a structure and 
the data essentially scatter symmetrically around the line Ip = I a which 
we display as a dashed straight line. This indicates that for the torus we do 
not find a pronounced pattern which links the localization of the mode to 
particular values of (. 

7 Summary 

In this article we have studied properties of topological excitations of SU(3) 
lattice gauge configurations with topological charge Q = ±1. Through the 
index theorem such configurations give rise to a single zero-mode of the 
Dirac operator. The zero-mode is localized on the topological excitation and 
reflects the properties of the underlying lump in the gauge field. We analyzed 
how the zero-mode changes when applying an arbitrary phase exp(i27r£) at 
the temporal boundary condition for the Dirac operator. Our main findings 
are: 

• For the ensemble in the deconfined phase we find that the zero-mode 
can change its position at £ = for configurations with real Polyakov 
loop and at Q = 1/3, Q = 2/3 for configurations in the complex sector. 

• The zero-mode is very spread out at these critical values and most 
localized for Q = 0.5 for real Polyakov loop, respectively C = 5/6, C, = 
1/6 for the complex sectors. 

• These observations nicely match the predictions for zero-modes of KvB 
solutions. 

• An analysis of the spectral gap shows the same periodic behavior in £ 
supporting earlier findings of only a single transition temperature for 
chiral symmetry restoration in all sectors of the Polyakov loop. 

• For the ensemble in the confined phase just below T c we find that for 
a large portion of our ensemble the zero-mode is located at different 
positions when changing the boundary condition. 
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• Also the localization of the mode fluctuates considerably during a com- 
plete cycle through (. 

• When analyzing the time dependence one finds that the zero-modes 
are stretched along the time direction for both ensembles above and 
below T c . 

• Several properties of zero-modes in the confined phase match predic- 
tions for KvB zero modes, but the resemblance is not as close as in 
the deconfined phase. 

• Also the zero-modes for configurations on the torus change their loca- 
tion when changing the phase at the boundary. 

• The results obtained when applying the boundary condition at differ- 
ent directions are consistent with each other. 

• About half of our torus configurations have zero-modes located at 
different positions when comparing periodic to anti-periodic boundary 
conditions. 

Our results show that for a large portion of our configurations an excita- 
tion with topological charge \Q\ = 1 is not a single lump. Instead it is built 
from several independent objects and does not resemble a simply connected 
instanton-like object. 
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